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ABSTRACT  Based  on  the  measured  altitude  distributional  data 
of  the  turbulence  structure  constant  (  C\  ) ,  this  article 
analyzes  the  effects  of  turbulent  outer  scale  Lo  on  atmospheric 
transverse  coherence  length  ro  and  isoplanatic  angle  0o  under  the 
two  conditions  of  large  apertures  and  small  apertures. 

Calculation  results  clearly  show  that,  under  conditions  where  Lo 
is  a  limited  value,  ro  and  ©o  clearly  increase.  With  regard  to 
large  aperture  telescopes,  in  infrared  wave  bands--for  example, 
wave  lengths  of  2.2  microns — ro  and  0o  are  even  larger,  ©o  can 
even  reach  0.5  mrad. 

KEY  WORDS  Turbulent  outer  scale  Atmospheric 
transverse  coherence  length  Isoplanatic  angle 
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1  INTRODUCTION 

Atmospheric  coherence  length  ro  and  isoplanatic  angle  0o  are 
two  important  parameters  describing  spacial  coherence 
characteristics  associated  with  light  wave  propogation  in 
turbulent  atmospheres.  Beginning  in  1968,  after  Fried  [1]  put 
forward  the  concept  of  coherence  length,  it  was  universally 
accepted  by  academic  circles.  In  recent  years,  research 
associated  with  ro  and  0o  has  grown  more  dynamic.  It  has 
already  spread  from  the  realm  of  laser  transmission  to  the  world 
of  astronomy.  The  source  of  the  impetus  lies  in  broad 
applications  associated  with  adaptive  optical  technology  (AOT) . 

As  is  widely  known,  AOT  is  capable  of  compensating  for  coherence 
characteristic  degnerative  regression  associated  with  light  waves 
when  they  pass  through  the  atmosphere.  Moreover,  the 
compensation  effects  are  directly  related  to  ro  and  ©o  .  ,  As  a 
result,  theory  developing  ro  and  0o  as  well  as  experimental 
research  are  very  important. 

It  is  generally  thought  that,  in  astronomical  observations 
and  measurements,  within  the  range  of  visible  light — ro  =  2-2 Ocm 
and  ©o  being  5-20microrad — when  telescope  apertures  D  are  smaller 
than  ro  ,  the  resolution  of  star  images  is  limited  by  the 
telescope.  When  D  is  larger  than  ro  ,  resolutions  are  limited  by 
the  atmosphere.  In  recent  years,  following  along  with  theory  and 
experimental  work  associated  with  in  depth  development,  there  are 
different  ways  of  looking  at  this.  This  is  primarily  the 
discovery  of  atmospheric  turbulence  scale  Lo  above  boundary 
layers.  This  certainly  does  not  seem  like  the  assumed  several 
tens  to  several  thousand  meter  ranges  in  traditional  theories, 
but  are  relatively  smaller.  If  Lo  is  actually  relatively  small, 
large  aperture  telescope  resolutions  will  be  greatly  improved. 

In  infrared  wave  bands,  it  is  even  possible  to  reach 
0 . 05micorrad.  At  the  present  time,  experiments  have  already 
clearly  shown  that — opting  for  the  use  of  wave  lengths  of  2.2 
microns — it  is  possible  to  reach  resolutions  which  are 
0 .25microrads  [2].  Recently,  Coulman  and  others  [3,4]  made  use 
of  balloon  borne  high  altitude  sounding  instruments  and  star 
light  scintillation  technology.  In  France,  the  U.S.,  and  Chile 
carried  out  wide  ranging  measurements,  guaging  Lo  values  to  be  30 
-  40cm.  Although,  in  the  vicinity  of  the  top  of  the  troposphere, 
Lo  was  measured  as  approximately  5m,  these  amounts,  however,  are 
not  adequate  to  influence  averages  associated  with  the  atmosphere 
as  a  whole.  McKechnie  [2,5]  believes  that  atmospheric  trubulence 
is  primarily  located  within  altitudes  4-5km  from  the  ground.  Lo 
associated  with  the  atmosphere  as  a  whole  are  only  30-40cm. 
Atmospheric  sounding  data  taken  from  Kunming  clearly  shows  that 
Lo  agrees  with  the  measurements  of  Coulman  and  others  in  order  of 
magnitude.  As  a  result,  it  is  then  necessary  to  have  the  effects 
of  finite  turbulent  outer  scales  on  atmospheric  coherence  lengths 
and  isoplanatic  angles  to  carry  out  research. 
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In  traditional  theory,  when  calculating  wave  structure 
functions,  one  generally  opts  for  the  use  of  Kolmogorov 
turbulence  spectra.  In  this  type  of  model,  it  is  believed  that 
Lo  is  guite  large,  with  the  result  that  it  may  not  be  considered. 
Now,  through  the  measurements  of  Coulman  and  others  as  well  as 
our  measurements  and  calculations,  it  was  discovered  that  Lo  is 
relatively  small.  The  equivalent  Lo  for  the  atmosphere  as  a 
whole  is  less  than  lm.  As  a  result,  when  calculating  wave 
structure  functions,  it  is  necessary  to  make  some  corrections  to 
Kolmogorov  spectra.  If  one  opts  for  the  use  of  deformation  Von 
Karman  spectra 

On  (K )  =0.033  C\  (K2  +  K2)~">*ex  p  (-  ~^T  ) 


Here,  K  is  the  spacial  wave  number.  The  unit  is  m-1  .  Ko  = 
1/Lo.  Km  =  5.92/lo  .  lo  is  turbulent  flow  inner  scale.  Cl 
is  a  characteristic  turbulent  flow  strength  index  of  refraction 
structural  constant.  The  unit  is  m— 2/3  .  Then,  wave  structure 
constants  [7]  are  „  „ 


D  (p,  z)  =  8 7r 2 k2 

J  0 


[W0  (K,p)]®,  (K)KdK 


(2) 


In  the  equations,  $n  (K)  must  make  use  of  equation  (1)  for 
substitution.  Moreover,  k=2nr/A.  X  is  wave  length.  Jo  (x)  is  a 
zero  order  Bessel  function.  p  is  the  distance  between  the  two 
points  considered.  In  order  to  make  calculations  convenient, 
take  $n  (K)  and  write  it  as  $n  (K)  =  $no  (K)  -  $1  (K,  Lo)  -  $2 
(K,  Lo,  lo).  Here,  there  is  no  relationship  between  $no  (K)  and 
Lo  or  lo  ,  that  is,  Kolmogorov  spectra.  $1  stands  for  the 
influence  of  Lo  . 

$2  stands  for  the  influences  associated  with  lo  . 


<3)  <D„0  (K )  =0.033  C2„  K~u,i 

(4)  <t>,  (K,  L„  )  =  0.033  C2n  (K2  +  K2)-"'6] 

(5)  (D2  (K,  L0,  /„)  =  0.033  C2„K-,,/3[l -exp  (-  )] 

Km 

In  equation  (6),  the  influences  of  Lo 
on  $2  are  ignored.  Going  through  detailed  analysis,  it  is 
possible  to  finally  obtain  [6]  p  <  lo  .  One  has 

r 00 

(6)  D  {p,  z)  =  3.28 /c 2 p 2  j  -gfrdz 

Jo  *  0 

For  lo  <  p  <  Lo  ,  one  has 


(7)  2.91  k2p512 


oo 

Cll  1  “■  0.8048  (  )in]dz 
Lo 
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For  Lo  <  p,  one  has 


relatively  large, 
spectra  to  derive 


1.563  k‘ 


(8) 


C\L?dz 


In  equation  (8),  when  Lo  is 
(p/Lo)  ->  0.  It  is  possible  to  use  Kolmogorov 


D(p,  z)=  2.91  k2p 5/3  C\dz. 


r 


3  INFLUENCES  OF  OUTER  SCALE  Lo  ON  ro  AND  ©o 

Use  was  made  of  balloon  borne 

atmospheric  sounding  data  taken  at  Kunming  by  Zeng  Zongyong  of 
the  Academia  Sinica  Anhui  Photomechanical  Institute  and  others  in 
order  to  analyze  the  influences  of  Lo  on  ro  and  ©o  .  Below, 
considerations  are  divided  between  two  cases.  One  case  is  lo  <  p 
<  Lo  .  The  other  is  p  >  Lo.  Through  another  type  of  expression 
for  wave  structure 

functions  [1]  /471 

D(p)  =  6M  (-*-  )5'3 

r«  (9) 


Making  equations  ( 7 )  and  ( 8 )  respectively  equal  to  equation  ( 9 ) , 
it  is  then  possible  to  obtain  (when  lo  <  ©  <  Lo  and  p  >  Lo) 
changes  in  ro  following  along  with  Lo  .  Then,  from  the 
relationship  [8]  ©o  =  0.314ro  /Ho  it  is  then  possible  to  solve 
for  corresponding  ©o  values.  At  this  time,  atmospheric 
equivalent  altitudes  Hq=[  Jcj  (z)zvij2/JC2  .  {z)dz]v^ 

Below,  respective  consideration  is  given  to  these 
two  types  of  cases. 

(1)  lo  <  p  <  Lo  .  Take  p  =  0.12m.  From 
equations  (7)  and  (9),  one  gets 

(10)  r0={  0.423/c2j*  C\ [  1  - 0.8048  (  -  ^ -  )1/3]<M~3/5 

The  altitude  distribution  associated  with  Lo  is 
taken  to  be  a  type  of  charted  formula  (see  Fig.l)  based  on  data 
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measured  in  Chile  [3] 


(11) 


L0  (z )  =  0.5  +  5  exp  ( -  Z  )2 


In  the  equation,  Lo  and  z  units  are,  respectively, 
meters  and  thouands  of  meters . 


Fig.  1 


Fig.  1  The  height  distribution  of  turbulent 
outer  scale 


To  make  comparisons,  using  the  formulae  derived  in 
association  with  Kolmogorov  spectra 

CD 

ro  =  [ 0.423  k2|  iz)dz]  3,5 


0O  =  [  2.9 1 A 2 


/•  CD 

•  C2n  (z)z5ndz]-}>5- 

Jo 


to  make  calculations,  the  calculation  results  for  the  two 
instances  are  seen  in  Tables  1A  and  2 . 

(2)  p  >  Lo  .  Take  p  =  8m.  From  equations  (8)  and  (9), 
one  gets  the  equation  below  (for  calculation  results,  see  Table 
IB)  . 


r0  =  {  0.227  A 2 


f  00 

C2„(z)(L0(z)]  ) 5/3 ] dz  } -3/5 p 
J  0 


(12) 
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TABLE  1 


Table  1  The  calculated  values  of  r0  and  0O  for 
Von  Kami  an  spectrum 


A.  p  =  0.12m  B.p  =  8m 


r0/ cm 

0o/prad 

'o 

/cm 

0O  /  prad 

x/^m 
time  X 

2.2 

0.55 

2.2 

2.2 

11  —  28  —  09 

7.8 

37.1 

2.18 

11.49 

30.7 

162.1 

9.51 

50.18 

11-28-15 

5.6 

29.4 

2.56 

13.54 

24.7 

130.6 

11.39 

60.13 

11-29-07 

20.0 

105.5 

8.64 

45.63 

166.7 

879.7 

72.09 

380.50 

11-29-19 

24.1 

127.4 

10.78 

56.91 

149.0 

786.5 

66.56 

351.30 

11-29-23 

20.9 

110.5 

11.99 

63.29 

167.4 

883.5 

75.90 

506.16 

11-30-  16 

17.7 

93.3 

5.49 

28.99 

88.4 

466.8 

27.49 

145.10 

“II  —  28  —  09*  means  at  9am,  on  Nov  .  28,  1992. 


TABLE  2 

Table  2  The  calculated  values  of 
Kolmogorov  spectrum 


p  =  0.12m 


'o 

/cm 

BJfiiad 

0.55 

2.2 

0.55 

2.2 

5.2 

27.6 

1.62 

8.54 

4.4 

23.2 

2.02 

10.67 

13.8 

72.6 

5.96 

31.45 

17.3 

91.2 

7.73 

40.81 

14.5 

76.4 

8.30 

43.79 

12.9 

68.0 

4.01 

21.17 

JxS'J 
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Considering  the  influences  associated  with  turbulent  outer 
scale  Lo,  use  was  made  of  expressions  for  ro  and  ©o  derived  with 
Von  Karman  spectra.  Making  use  of  empirically  measured  C\ 
vertical  distribution  data,  calculations  were  carried  out  with 
regard  to  ro  and  ©o  .  From  the  calculation  results,  it  was 
discovered  that,  after  considering  the  influences  of  Lo  ,  ro  and 
0o  clearly  increase.  In  the  past,  it  was  generally  believed  that 
there  was  practically  no  difference  between  the  performance  of 
large 

aperture  telescopes  due  to  atmospheric  turbulence  and  small 
aperture  telescopes  of  10-20cm.  Now,  considering  Lo  effects, 
estimates  are  much  larger  making  use  of  large  aperture  telescope 
ro  and  ©o  than  using  traditional  theories.  In  particular,  in  the 
infrared  wave  length — for  a  2.2  micron  aperture,  for  instance — it 
is  8m,  and  0o  is  capable  of  reaching  approximately  0.5mrad.  With 
regard  to  tracking  and  guidance,  this  is  capable  of  having  a  much 
larger  selection  domain.  ro  is  capable  of  increasing  to  being 
equivalent  to  apertures  in  contact.  In  this  way,  the  receiving 
system  is  capable  of  possessing  optimum  properties  [9].  In  this 
type  of  case,  adaptive  optical  technology  is,  perhaps, 
unnecessary,  or  adaptive  optical  systems  with  relatively  low 
indices  are  also  capable  of  satisfying  requirements.  We  expect 
that  this  conclusion  can  be  experimentally  verified. 
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ABSTRACT  On  the  basis  of  transient  state  four  dimensional 
(4D)  thermal  blooming  programs,  this  article  studies  transient 
thermal  blooming  effects  during  atmospheric  laser  transmission  as 
well  as  how  they  transition  to  thermal  blooming  states  under 
isobar  approximation.  In  conjunction  with  this,  comparisons  are 
made  with  4D  thermal  blooming  program  calculation  results  under 
isobar  approximation.  At  approximately  four  fold  fluid  dynamic 
time  scale  (4th),  transient  4D  thermal  blooming  program 
calculation  results  are  close  to  results  under  isobar 
approximation.  Besides  that,  phase  compensation  calculations 
were  also  done.  In  cases  with  light  intensities  of  a  few  kW/cm2, 
under  ideal  adaptive  optical  system  kinetic  conditions,  transient 
phase  thermal  distortions  are  capable  of  being  compensated  for  in 
real  time  relatively  effectively. 

KEY  WORDS  Transient  thermal  blooming  Isobar  approximation 
Phase  compensation 
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I .  BASIC  EQUATIONS 

Thermal  blooming  effects  are  limitations  on  laser  strengths 
passing  through  the  atmosphere,  giving  rise  to  nonlinear 
phenomena  associated  with  beam  distortions.  Due  to  the  fact  that 
the  details  of  theoretical  research  outside  China  on  thermal 
blooming  effects  during  the  full  course  of  laser  transmission  in 
the  atmosphere  have  still  not  been  made  public,  as  a  result, 
there  is  a  great  need  for  research  on  laser  thermal  blooming 
effects  during  the  whole  process  in  the  atmosphere  as  well  as  for 
the  establishment  of  4D  thermal  blooming  programs .  This  article 
lays  stress  on  the  description  of  transient  effects  in  order  to 
facilitate  an  understanding  of  whether  or  not  transient 
distortions  are  capable  of  effective  real  time  compensation.  In 
conjunction  with  this,  observations  are  done  of  the  effectiveness 
of  isobar  approximations.  Use  is  made  of  computer  numerical 
values  to  describe  basic  equations  associated  with  thermal 
blooming  of  the  several  types  below. 

1.1  Field  Equations 

Using  monochromatic,  linear  polarized  scalar  wave  motion 
equations  to  describe  laser  fields — under  asymptotic 
approximation — in  following  the  wave  motion  coordinate  system, 
the  slowly  changing  field  propogated  in  the  positive  Z  direction 
satisfies  the  equation  below 

2ik-^ +V2±ip  +  k25e(il/,  t)\j/  =  0  (1) 

In  this,  k=2rc/A.  A  is  laser  wave  length  in  a  vaccuum.  is 

the  Laplace  operator  associated  with  vertical  light  axses.  6e 
is  diaelectric  constant  change. 

As  far  as  beacon  light  propogated  in  the  negative  Z 
direction  is  concerned,  it  has  the  equation  below 
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-2ikb-^-  +  Yi  (2) 


Assuming  that  kb  =  k  ,  6eb  =  6e(i|i,t),  and  omiting  delay 
effects,  it  is  believed  that  beacon  light  is  reverse  propogated 
in  a  "distribution  lens"  created  by  laser  emissions.  Nonlinear 
Xue  Dingwo  equation  (1)  can  use  various  types  of  methods  for 
solution  (integral  transform  methods,  difference  methods,  as  well 
as  Green  Function  methods,  and  so  on).  In  the  processing  below 
(with  a  0( AZ2 )  precision)  using  FFT  methods  to  solve,  one  takes 
equation  ( 1 )  and  solves  in  terms  of  a  form  within  a  step  length 
of  AZ  to  get 

V*  ItSF 

HZ  +  AZ,  r_i ,  f)  =  exp[i(  )AZ]\j/(Z,  r±  .  r)  (3) 


1 

dE=  - 


Se(  Z')dZ' 


r\ 


In  this 


AZ 


1 


is  the  vertical 


optical  axis  coordinate  vector. 


Equation  (3)  can  be  expressed  as 

iMZ  +  AZ,  r  ,  r)  =  exp[i(-^-)  -y-  ]exp[i(— —  )AZexp[i(-^-  )-^-]*KZ,r±,  t) 

When  making  use  of  equation  ( 3 '  ) ,  it  requires  step  length 
AZ  to  satisfy  the  condition  below  [1] 


(4) 


~  SsAZ«  1 
2 


far  as  using  the  following  structural  relationships 
as  well  as  density  changes  in  isobar  approximations  are 
concerned,  with  regard  to  propogation  calculations  for  laser 
beams  with  wave  lengths  of  1  micron,  it  is  necessary  for  step 
length  AZ  to 

satisfy  , 

AZ/cm«  — — -  =2.5  x  10  8 — — TTTT^Z^WtJl  (4') 

\p,\  a  (cm  ')/./(MWcm  )'t/ s_  v**  ) 


In  this  relationship,  aa  is  the  absorption  coefficient 
peak  value  light  strength,  t  is  time. 

~  _ 1 _ 1 _  T7 


Ip  is 


Obviously, 


V?  x  AZ  .  , 
exp[i(  -rf-  )  ]^(2,r±,  t) 


in  equation 


(3-)  is  a  solution  in  the  back  half  of  step  lengths  for  Z  in  the 
vaccuum  wave  motion  equation  below 


2  ik 


d\f/ 

~dZ 


+  Vy  tp  =  0 


(5) 


The  effects  of  dielectric  media  are  equivalent  to  phase^ 
multiplication  factor  exp[i(y)^AZ],  *  In  t  lS 

fields  obtained,  use  is  then  made  of  operator  exp^-2T')“2-^ 

to  act  on  them.  This  is  equivalent  to  taking  the  obtained  fields 


(in 


7+  locations)  and  making  "initial "  conditions  on 

2 

the  bottom  half  of  step  lengths  to  solve  vaccuum  wave 
equations.  In  this  way,  one  gets  a  step  length  field.  In  the 
same  way,  step  lengths  are  solved  for,  one  after  the  other 
(called  multiple  phase  screen  treatment),  that  is,  it  is  possib 
to  solve  for  fields  at  different  propogation  distances. 

As  far  as  solving  for  vaccuum  wave  motion  equations  in 
configurational  space  is  concerned,  in  solutions,  there  are 


X 
/  ^ 


) 


cP 

contained  second  degree  differential  operators  V”x(  --&+T? 
Moreover, 

in  " spaciai^f  requency  space,  then,  from  eguatron  (5),  one  has 

^(Z  +  AZ,  fcx,fcy)  =  exp[-  ^-(^  +  ^)AZ]^(Z,fex,fcj,)  (6) 


Here, 

is, 


>/> 


is  the  Fourier  transform  associated  with 


/•A*00 


i|r  ,  that 
/375 


MZ.kx,ky)  = 


i//(Z,  rL)exp(-iA*r1  )^r± 


Besides  this. 


= jj  i/f(z,  x,  y)exp[-i(ferx'+fc),y)]^>:<iy 
r  r°° 

\f/(Z,X,Y)  =  -^T  I  ${Z,kx,ky)e*v(ik-r±)dkxdky 


conjunctio^with  boundary 

conditions  to  obtain 


,  _  2n  i  -  2  71  - 

fc*  m  ,  ky  .  n 

Lx  *- y 


-  W/9  12  N/2  and  Lx,  Ly  are  function  periods 

In  this,  m,n--N/2. . *1'2'**  (  the  two  to  be  equal,  it  is 

in  the  two  directions.  Choosing  the  twc^  transverse 

necessary  to  point  out  that,  if  are  steep,  then,  there 

distributions  in  configurable Contained  in  "special  frequency" 
should  de  a  great  ^any  large.  and 

configurationalS space  coordinate  disperse  shows 

X<=iAX  ,  Y,  =  IAY 


.  ,  _  i  M  -i  Takina  transverse  step  lengths 

In  this,  D / 1-0 JJ  h  made  of  the  various  equations  above. 
AX=AY=L/N,  application  is  maae 

Then,  equation  (6)  becomes 

^m„(Z  +  AZ)  =  exp[-  -±-{^-)2(m-  +  n2)AZ]iPmn(Z)  (6') 

Thus ,  optical  strengths 


/(Z  +  AZ,  Xj,  7,)=  ~^M\Z  +  AZ,  Xj,  Y,) 


ti  \ 


Giving  consideration  to  linear  absorption  and  scattering, 
there  is  a  need  to  introduce  attenuation  operators 


exp[- 


/•Z  +  AZ 

*t(Z')dZ'  ], 


JZ 


at  -  as  +  aa  .  aa  and  as  are,  respectively,  absorption  and 
scattering  coefficients. 

1.2  Fluid  Dynamic  Equation  Sets 

Due  to  the  fact  that  atmospheric  laser  absorption  for  wave 
lengths  that  interest  us  is  very  weak,  during  laser  transmission, 
it  leads  to  atmospheric  status  parameter  changes  being  very 
small.  It  is  possible  to  opt  for  the  use  of  linearized  fluid 
dynamic  equation  sets  to  describe  changes  in  atmospheric 
configuration  parameters,  that  is. 


dp, 

+PoV*  u  =  0 

(8a) 

du 

Po - =  -cjVp, 

(8b) 

dt 

2  d£i  _  ■>  dp,  _ 

2  dt  3  dt  ^ r 

(8c) 

In  these,  because  pressure  pi  units  make  use  of 

atm( latm=101 . 325kPa)  and  optical  strength  Ip  units  make  use  of 

MW.cm-2  ,  and,  due  to  introduced  unit  transformation  operators 

q=  1 .0133  x  106,  c;  =  0. 10133  x  10-6,  <~=  10_l3c;,  ,  r  is  an 

adiabatic  index,  and  Ip  does  not  depend  on  time  t.  It  is 
possible  to  obtain  scattering  spectra  associated  with  density 
changes  and  pressure  changes  in  "spacial  frequency"  space 

PUZ,t)=-  -M^-rfl-  jjp(W^2  +  «Mi  (9 


Plmn(Z  )f)~ 


Sm„(Z) 


nr  +  rr 
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■sin(wtyfm*+n*  )  (10) 


In  this,  \v=  — 1  —  =c  2a  ,  and 

c2  L  s  L 

Sm«(Z)  =  (r-  l)aaexp[- 

Doing  reverse  Fourier  transforms,  it  is  possible  to  obtain  pi  and 
pi  on  scattering  points  in  configurational  space 


a  ,{Z')dZ'}lpmn(Z) 


(11) 


Pi(Z,j Ax,  lAy ,  t)  = 


X  X  P\mn(Z ,  f)exp[i  -^j-(m/+n/)J 

,n=  ~N/2  iV 


(12) 


(13) 


n/2  -I  ^ 

Pi(Z,jAx,  lAy,  t)=  £  £  plmn(Z,  t)exp[i  —rr(mj+nl)] 

m.n=-A72  iV 

Applying  ideal  gas  configuration  equations,  it  is  possible  to 
obtain  temperature  changes 

T\(Z ,jAX,  IAY ,  t)=  -Q—  —  T0—  (14) 

PoF  Po  1  ' 


In  this,  r=2 . 848  atm. cm3/K.g. 

From  equation  8(c),  it  is  known  that  isobar  conditions 
correspond  to  c2  =  0.  When  Ip  does  not  depend  on  time  t, 
integrating  t,  one  gets  density  changes  associated  with  isobar 
conditions 

i>]p(Z,X,Y,t)=-±S(Z,X,Y)  . 

ci  (15) 

This  and  the  results  in  reference  [2]  are  the  same.  Making 
Fourier  transforms  with  regard  to  equation  (15),  scattering  gives 
one 

plmnp(Z,X,Y,t)=-±Smn(Z)  (16) 


This  is  precisely  the  first  term  on  the  right  side  of 
equation  (9).  From  this,  it  is  possible  to  know  that  transient 
correction  terms  at  nonisobaric  times  are  precisely  the  second 
term  contained  in  square  brackets  on  the  right  side  of  equation 
(9).  They  follow  increases  in  t  and  decrease.  Moreover,  high 
frequency  mode  contributions  follow  increases  in  t  and  attenuate 
even  faster.  With  regard  to  Gaussian  distributions,  low 
frequency  modes  account  for  the  dominent  positions.  In  this  way, 
when  Xm  =  10,  from  sinX/X  ->  0,  it  is  then  possible  to  estimate 
that  transient  contribution  attenuation  times  are  a  few  fluid 
dynamic  time  periods  (Here,  fluid  dynamic  time  scales  tH  =  a/cs. 
a  is  optical  strength  e-2  beam  radii.  Takinq  transverse  function 
period  L=6a,  attenuation  time  period  td=Xmt, t/j  m2+n2.  )  • 

Besides  that,  from  equations ' containing  source  souna  waves,  it  is 
also  possible  to  derive  this  transient  term.  From  (8a  -  8c),  it 
is  known  that, 


1  <52  _ 

c:  or  r‘ 


r-1 


(17) 


Isobar  approximations  correspond  to  the  speed  of  sound  cs  -> 
oo,  that  is,  distrubances  in  the  development  of  infinite  speeds, 
instantaneously  bringing  them  to  the  same  level,  when  Ip  does  not 
depend  on  time  t.  From  equation  '27'*  it  is  possible  to  obtain — 
under  isobar — density  Pip=  ~ -ip— “aV*  •  This  is  the  same 

as  equation  (15).  In  early  period^,  pi  transverse  changes  are 
small.  V|p,-*  0,  •  When  Ip  does  not  depend  on  time  t,  it 

is  possible,  from  equation  (17),  to  obtain  density  variations 


(18) 


associated  with  t3  thermal  blooming  phases 


Normally,  wind  speed  vw  «  cs  .  Thus,  wind  transit  periods 
are  far  longer  than  fluid  dynamic  time  scales  tH  .  As  a  result, 
thermal  blooming  phenomena  within  fluid  dynamic  time  period 
scales  for  research  are  capable  of  not  figuring  in  the  influences 
of  wind.  In  this  way,  equation  (17)  is  axially  symmetrical.  /377 
When  Ip  does  not  depend  on  time  t,  it  is  possible  to  go  through 
Han  Keer  transforms,  and  that  will  do.  With  regard  to  second 
order  ordinary  differential  equations  associated  with  °pxfit  , 

use  is  made  of  initial  conditions  ,(_q  dPl 

/r_n  \  =  q  .  It  is  possible  to  obtain^'  ’  -  ~  ~W 


Pi(t,r±,  z)  = 


Ip(t)sm(csXt)*  J0(Ar±  )d'A  19) 
o 


In  early  periods,  it  is  thought  that  Ip  does  not  underao  thermal 
blooming  effects.  It  is  possible  to  adopt  Ip(r^_,  :)  =  I°pex p 

(-2| r  1 2/a2),  •  Making  Han  Keer  transformations,  one  obtains 

/PU)=  -fa2exp  ~) 

On  optical  axses,  r±  =  O.  Thus,  zero  order  Bessel  function  Jo 
=  1.  Because  of  this,  at  a  relatively  early  time,  from  equation 
(19),  one  obtains 

pl(t,r±  =  0.Z)  =  -  -^r-aa^(r-  -±L)  (20) 


From  this,  it  can  be  seen  that  a  few  fluid  dynamic  time  period 
transient  portions  are  then  attenuated  away. 

1.3  Structural  Relationships 

Structural  relationships  give  the  relationships  between 
indices  of  refraction  or  dielectric  constants  and  dielectric 
media  configuration  parameters.  From  reference  [2],  one  knows 
that 


«-  1  ~  —  X  10 —  p  =  0.23p 
Po 


(21) 


From  this,  it  is  possible  to  obtain 


<5e=  2dn  =  0.46/9, (Z,  A",  Y,t) 


2  CALCULATION  RESULTS  AND  ANALYSES 


(1)  Definite  Solution  Conditions  (Collimated  beam, 

Gaussian  transverse  distribution) 


Z  =  0,  <j,(X,Y,Z  =  0)=^0exp(-^!±il) 

cr 


(22) 


In  this, 


iAo=2.71  xio4/^ 


.  If  consideration  is  given  to 


phase  compensation,  the  right  side  multiplication  factor 


,exp{-itg-'[ 


ImiM*,  r.z-o)  n 

Re<l>b(X,Y,Z=0)  J' 


,  z> o  jr.y— 


±cc, 


0,  Z  —  Lg  ,  MX,  Y,Z=L„  )  =  ^0exp(-  ) 


In  this,  A/  is  laser  emission  transmission  distance,  t  =  o. 
pl=pl=u=0.  (2)  Atmospheric  Model  p(Z)oce-«- 

b  =  0 . 143xl0-5cml .  (3)  Parameters  Used  A  =  1  micron, 

aa(Z=O)=0. 12xl0-6cm-l,  at(Z=0)=0. 12xl0-6cm-l,  r=1.4;  cs  =  330m/s. 

Use  was  made  of  transient  thermal  blooming  4D  programs 
(FFT  methods)  to  do  numerical  value  calculations  on  thermal 
blooming  associated  with  different  beam  radii  ao  ,  different 
optical  strengths  fp,  ,  and  different  transmission  distances 
.  Results  are  seen  in  the  various  tables  or  charts.  For  the 
sake  of  comparison,  numerical  value  calculations  associated  with 
isobar  approximations  were  done  on  certain  models  (Table  1).  In 
conjunction  with  this,  phase  compensation  calculations  were  made 
(Tables  3,  4). 

From  Fig.l,  it  can  be  seen  that,  when  t  <  2tH,  pressure 
change  radial  distributions  are  basically  pulled  to  the  same 
level.  After  4tH,  an  even  status  was  maintained  throughout. 
However,  beam  edge  areas,  by  contrast,  are  higher  than  beam 
centers.  This  is  due  to  the  use  of  continuous  wave  Gaussian 
radial  distribution  lasers.  In  this  way,  center  location 
pressure  perturbation  sources  are  maintained  right  along. 

Pressure  disturbance  waves  are  propogated  toward  beam  edges.  Due 
to  "pursuing  characteristic  line"  effects  (characteristic  line 
converegence ) ,  edge  pressures  are  slightly  higher  than  the 
center.  Besides  that,  it  is  also  possible  to  see  that,  with 
regard  to  beam  center  vicinity  fixed  spacial  positions,  pressure 
perturbation  value  time  changes  go  up  first.  They  do  not  reach 
tH  and  then  begin  to  decline.  When  they  reach  the  vicinity  of 
4tH,  they  do  not  decline  any  more.  After  that,  they  gradually  go 
up  again.  The  explanation  of  this  is  constant  source 
perturbation.  With  regard  to  time  scales  being  unable  to  form 
real  isobar  configurations,  although  pressure  space  distributions 
are  capable  of  being  pulled  up  to  the  same  level,  fixed  point 
pressures,  however,  still  follow  along  with  increases  in  t  and  go 
up.  This  type  of  so  called  "isobar  configuration"  is,  in  fact,  a 
superposition  of  a  ^umber  of  transient  effects.  Giving 
consideration  to  #o,  ,  pi  will  be  smaller  than  equation 

(15)  when  at  times  ui-'loubaric  approximation  (see  Table  1  last 
column) .  /378 
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TABLE  1 


m  1  fl=250an,  J,-3.2kW  •  cmT2 

Table  1  The  on  axis  parameters  obtained  by  various  methods 


L,,  =  3km 

5km 

7km 

10km 

r/ms 

h 

~A 

Pi 

h 

“A 

P\ 

h 

-A 

P\ 

h 

-A 

Pi 

A 

2.60 

1.100 

2.30 

0.710 

2.10 

0.490 

1.90 

0.290 

1.5 

(0.2 tH) 

B 

2.60 

0.056 

1.1 

2.3 

0.038 

0.76 

11 

0.026 

0.53 

1.90 

0.016 

0.016 

0.32 

C 

0.058 

0.039 

0.027 

3.75 

A 

2.55 

2.600 

2.29 

1.80 

2.10 

1.20 

1.90 

0.72 

B 

2.55 

0.740 

2.1 

2.30 

0.50 

1.50 

110 

0.35 

LQ 

100 

0.20 

0.25 

0.60 

(0.5 t„) 

C(18*) 

0.900 

0.61 

7.5 

A 

2.55 

5.20 

2.30 

3.50 

2.00 

140 

1.80 

1.40 

B 

2.55 

3.80 

2.00 

2.30 

2.40 

1.30 

100 

1.70 

0.90 

1.90 

0.95 

0.55 

(U») 

C 

15 

A 

2.50 

10.0 

| 

2.20 

6.70 

100 

4.60 

1.80 

170 

B 

2.50 

10.0 

0.80 

2.20 

6.50 

0.56 

100 

4.50 

0.39 

1.90 

160 

0.23 

(2  tH) 

C(15‘) 

12.0 

8.10 

5.60 

3.30 

30 

A 

2.49 

21.0 

2.00 

13.0 

1.90 

9.60 

1.60 

4.80 

B 

2.49 

21.0 

0.39 

2.00 

12.0 

0.26 

1.80 

8.40 

0.18 

1.50 

4.50 

0.11 

(*„) 

C(  15*) 

24.0 

16.0 

13.00 

6.80 

Note  :  Calculation  parameter  a=  250cm,  iJ)  =  3.2kW*  cm"2,  lp  is  in  kW-cm  2,  pressure  is 
in  10”6atm,  density  p,  is  in  10"9g  *  cm"3,  A  is  isobaric  results,  B  is  transient,  C  is 
analytical  solutions  (see  eq.  18),  *  Digit  after- C*  is  equation  number. 


From  Table  1,  it  is  possible  to  see  that,  at  different 
transmission  distances,  changes  in  (-pi)  as  a  function  of  t  are 
basically  in  agreement  with  transient  analysis  results  associated 
with  equation  (18)  when  t  <  tH/2.  When  t  >  2tH  ,  there  is, 
basically,  agreement  with  isobaric  approximation  results 
atlociated  with  equation  (IS).  When  the  two  really  approach  each 
r.+'hor'  ic  after  t  >  4tH.  On  axis  pressures  follow  changes  in  t 
not  with  monotomic  rises  to  saturation  values,  but  with  decreases 
after  first  rising,  and  with  continuous  waves  slowly  going  up. 
Extreme  value  time  periods  to,  can  be  calculated  from  equation 
(8c).  Adopting  -^-(r0)=°  '  thus' 

'  at 


lu{to)-~~k{r  w'° 


Taking 

(18). 


lp=Ip  exp(-2r|  jar) 


substitute  into  equation 


ri-o 

_  l)afl 

c\ar 


in  conjunction  with  that,  when  attention  is  paid  to 
n  '  8 /•  .  Thus,  4p,(t0)  = 

dt 


-  -77 

Because  of  this. 


tr\  — 


4Cj 


4, 


(23) 
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Moreover,  time  periods  when  on  axis  pressure  changes 
extreme  values  associated  with  transient  state  4D  thermal 
blooming  program  calculation  results  are  approximately  0 . 6  . 

The  ?wS  a?S  basically  in  agreement.  Rough  estimate  values  for  on 
pressure  changes  at  U  -  10km  locatrons  are 
— w  i  p, (r |  =  0,  =  1 0km,  to) —  P\ (r±~ 8 .  Z-0,  tn)  exp 

d .o oxiu- 1 atm® 0 . 5 68Pa  and  are  also  in  oasic  agreement  with  program 

calculation  results.  Here 


e-bzdz= 


(l-e  bL‘  )  (24) 


*«* 


g,(0) 

L‘i  Jo 


«t(0) 

bL„ 


Table  2 


Table  2  The  results  computed  by  transient  thermal  blooming  code  for  various  a  .  $  and  L.  at  4 1„ 


a/cm 

f 

Li,  /km 

b 

“ P\ 

Pi  F 

’=  kcr!L„ 

-V0 

3.2 

5.2 

0.27 

2.6 

1 

10 

32 

0.1 

16 

0.8 

8.1 

6.28 

47 

0.024 

26 

0.35 

320 

257 

0.0013 

266 

1.80 

0.4 

0.337 

0.66 

10.1 

1.91 

50 

3.2 

2 

2.43 

0.048 

81 

V~) 

CO 

13.9 

32 

20.4 

0.0040 

810 

Ill 

3 

2.49 

0.0021 

390 

13090 

91 

250 

3.2 

5 

2.0 

0.0012 

260 

7862 

7 

1.8 

0.084 

160 

5616 

10 

1.5 

0.045 

110 

3931 

Note  :  The  phase  compensation  have  been  unstabilty  when  ND  is  more 

than  several 

hundreds 

/•Lf 

then  N 

bndzjk  ,  Nc  = 

o 

=  Nd/F  (see  proc  SPIE  Vol.1221),  ij1 

,  Ip  are  in  kW  •  cm"2, 

~Pi  is 

in  10' 

'9atm,  p,  is  in 

io-9g  • 

cm  \ 
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From  Table  1,  it  is  also  possible  to  see  that,  with  regard 
to  fixed  instants,  pi  and  pi  follow  changes_in  transmission 
distances  L  „  and  are  basically  exp[- (6  +  a,)f 

Because  pi  and  pi  are  in  direct  proportion  to  alp,  of  course,  if 
thermal  distortion  is  severe,  Ip{Ltl)  cc  (-a,L„)  •  exp(-Nc) 

With  regard  to  our  case,  distortion  is  weak,  exp(-Nc)->l . 

From  Table  2,  one  knows  that,  no  matter  whether 
radii  are  large  or  small,  optical  strengths  are  high  or  low 
(within  a  certain  limit),  transmission  distances  are  long  or 
short,  thermal  distortion,  when  t  £  4tH,  is  very  weak  in  all 
cases.  Moreover,  on  axis  optical  strengths  and  h  form  a 


TABLE  3 


Table  3  The  calculation  results  for  phase  compensation  (radial  distribution  of  intensity) 


No 

65  (r  =0) 

66  ( 12cm) 

67  (24cm) 

68  (36cm) 

69  (36cm) 

70  (36cm) 

he 

1.78 

1.77 

1.75 

1.71 

1.65 

1.58 

1.97 

1.96 

1.94 

1.90 

1.84 

1.76  ■ 

h 

3.2 

3.19 

3.14 

3.08 

2.98 

2.87 

i; 

1.972 

1.966 

1.935 

1.898 

1.837 

1.769 

No 

71  (36cm  ) 

72  (36cm) 

73  (36cm) 

74  (36cm) 

75  (36cm) 

76  (132cm) 

he 

1.50 

1.41 

1.31 

1.21 

1.10 

1.00 

h'e 

1.68 

1.59 

1.44 

1.38 

1.27 

1.16 

h 

2.73 

2.58 

2.42 

2.24 

2.06 

1.88 

K 

1.682 

1.590 

1.490 

1.380 

1.269 

1.159 

Note  ;  “No*  is  radial  coordinate,  Ipc{I^)  is  intensity  without  (or  with)  compensation,  Ig(r^  )  is 
intensity  of  Gaussian  distribution,  /„  is  intensity  without  bloming,  all  of  intensity  units 
are  in  kW-cm-2,  r=4iH=30ms,  L„=10km,  a=250cm,  J,°=3.2kW- cm'2,  V  =  /g(r 
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direct  proportion.  This  is  extremely  beneficial  with  regard  to 
real  time  phase  compensation. 

In  Table  3,  thermal  distortions  associated  with  t  ^  4tH 
phases  are  not  great.  Compensation  is  easy.  Thermal  distortion 
optical  strength  radial  distributions  deviate  from  Gaussian 
distributions.  After  going  through  phase  compensation,  10km 
location  optical  strength  radial  distributions  almost  change  into 
ideal  Gaussian  distributions. 

From  Table  4,  it  is  possible  to  know  that,  no 
matter  how  early  the  time,  at  different  transmission  distances, 
very  good  compensation  was  obtained  in  all  cases.  Here  are  the 
calculations  carried  out  with  regard  to  continuous  lasers.  We 
have  not  yet  done  calculations  for  pulse  lasers .  Provided  pulses 
are  smaller  than  fluid  dynamic  time  scales — basically  belonging 
to  t3  thermal  blooming — it  is  possible  to  go  through  real  time 
phase  compensation,  correcting  transient  thermal  blooming 
distortion.  With  regard  to  pulse  lasers,  when  pulse  internals 
are  larger  than  or  egual  to  wind  transit  periods,  pulses  series 
will  not  accumulate  thermal  distortions  to  a  level  which  cannot 
be  compensated  for.  Reference  [3]  already  pointed  out  that 
single  pulse  and  multiple  pulse  thermal  blooming  are  both  weaker 
than  is  the  case  with  continuous  wave  lasers.  Therefore,  with 
regard  to  thermal  blooming  associated  with  pulse  series  laser 
beams,  it  seems  even  easier  to  make  adaptive  optical  compensation 
in  real  time. 
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Fig  1  The  radial  distributions  of  pressure  changes 
at  various  times  when  L„=100m 
(a=  1cm,  ;,°  =  32kW  -cnT2) 

3  PRELIMINARY  CONCLUSIONS 

In  order  to  save  on  computation  time— -time  when  computers 
are  in  use — spacial  step  lengths  are  relatively  large.  As  a 
result,  the  precision  of  calculation  results  is  not  too  high. 
However,  a  number  of  patterns  revealed  are  still  credible.  The 


reason  is  that  field  calculations  in  programs  go  through 
inspections  of  round  vaccuum  aperture  refractive  analysis 
results,  and,  in  the  programs,  the  pi  and  pi  obtained  from  fluid 
dynamic  calculations  agree  with  rough  analytical  estimate  values 
in  cases  of  limitations.  On  the  basis  of  the  calculations  above, 
it  is  possible  to  obtain  the  preliminary  conclusions  below: 

(1)  Transient  effects,  up  to  approximately  4tH  times, 
basically  approach  quasi  isobaric  states.  In  the  case  of  t  <  tH 

,  one  is  placed  in  a  t3  thermal  blooming  state.  In  the  case  of  t 
k  4tH  ,  one  basically  reaches  a  quasi  isobaric  situation.  When 
optical  strength  changes  are  not  considered,  density  changes  and 
time  t  form  a  direct  proportion. 

(2)  Within  the  range  of  parameters  which  interest  us,  no 

matter  how  large  or  small  a  is,  and,  no  matter  how  long  or  short 
1  is  or  how  high  or  low  optical  strengths  sre,  provided  one  is 
within  the  period  ,  thermal  distortions  are  all  very  weak. 

Thus,  it  is  possible  to  realize  relatively  effective  phase 
compensation.  However,  here,  the  effects  of  small  scale  /381 
instabilities  are  not  figured  in.  Moreover,  phase  compensation 
is  ideal. 

From  transient  4D  thermal  blooming  program  calculation 
results,  it  is  possible  to  predict  that,  with  regard  to  laser 
thermal  distortions  for  systems  associated  with  pulses  smaller 
than  fluid  dynamics  periods  tH  it  is  possible  to  do  real  time 
phase  compensation.  Provided  pulse  intervals  are  larger  than 
wind  transit  times — with  regard  to  requirements  associated  with 
"optimum"  characteristics  of  pulse  interval  structures — it  seems 
to  be  possible  to  relax  them.  Moreover,  in  regard  to  pulse 
spacial  structures,  there  is  still  a  need  for  attention  because, 
during  t3  thermal  blooming,  p,°cV%  .  As  far  as  concrete 

conclusions  are  concerned,  they  await  the  next  step  in  the 
detailed  analysis  of  4D  thermal  blooming  program  numerical  value 
calculations  concerning  such  effects  as  calculated  turbulent 
flow,  wind,  and  so  on. 
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